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INTRODUCTION 


This  is  a  second  annual  report  on  collaborative  study  between  AeroChem  Research 
Laboratories  (principal  investigator  -  Dr.  H.  F.  Calcote)  and  the  Pennsylvania  State  University 
(principal  investigator  -  Dr.  M.  Frenklach).  The  ultimate  objective  of  this  program  is  to  develop  a 
quantitative  physical/chemical  model  of  the  mechanism  of  soot  formation  that  is  consistent  with 
available  experimental  data.  The  specific  objectives  of  the  proposed  3-year  study  are:  1)  to 
delineate  the  relative  importance  of  radical  and  ionic  mechanisms  of  soot  nucleation  in  flames,  and 
2)  to  determine  the  optimum  chemical  and  physical  model  of  the  total  soot  formation  process.  Only 
the  Penn  State  work  is  reported  here,  that  of  AeroChem  IS  reported  in  their  annual  report1 

WORK  PROPOSED 

The  specific  objectives  proposed  for  the  second  your  of  the  of  the  Penn  State  efforts  were: 

1.  To  undertake  and  complete  the  analysis  of  the  flame  simulations  with  the  ion-reaction 
mechanism  provided  to  us  by  AeroChem. 

2 .  To  complete  the  overall  analysis  of  the  competition  of  the  PAH  formation  between  the  ionic  and 
radical  mechanisms. 

3 .  To  review  the  reactions  and  the  associated  rate  and  thermochemical  data  for  the  production  of 
larger  species  and  soot  particles  developed  by  AeroChem. 

4 .  To  develop  a  computer  code  for  modeling  the  formation  and  growth  of  large-sized  species  and 
soot  panicles. 

5 .  To  perform  computer  simulation  of  soot  particle  formation  in  flames  with  the  code  developed. 

RESULTS 

During  the  last  twelve-months  period  of  the  project,  we  completed  most  of  the  planned  work: 
the  ion  reaction  mechanism  provided  to  us  by  AeroChem  was  “cleaned-up”  and  run  with  a  flame 
code,  the  computational  results  were  analyzed,  a  computer  code  for  modeling  large  species  and 
soot  particle  formation  and  growth  was  developed,  our  neutral  reaction  mechanism  was  updated 
and  tested,  and  a  detailed  computer  simulation  of  soot  particle  nucleation  and  growth  in  laminar 
premixed  flames  was  successfully  performed  using  our  updated  neutral  reaction  mechanism.  The 
latter  part  is  a  particularly  significant  accomplishment,  as  it  is  for  the  first  time  soot  particle 
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inception  could  be  modeled  from  first  principles,  starting  with  fuel  decomposition  and  going  all  a 
way  to  predicting,  in  quantitative  agreement  with  experiment,  soot  particle  properties. 

Ion  Reaction  Mechanism 

The  flame  simulation  using  the  mechanism  provided  to  us  by  AeroChem2  was  completed.  The 
initial  results  indicated  that  a  number  of  the  rate  coefficients  are  greatly  overestimated  at  low 
temperatures,  causing  severe  numerical  problems  with  the  flame  simulations.  This  problem  was 
corrected  by  estimating  the  rate  coefficients  for  the  exothermic  direction  and  by  calculating  the  rate 
coefficients  for  the  reverse,  endothermic  direction  from  the  principal  of  detailed  balancing,  as  was 
advocated  in  our  work  on  the  neutral  reaction  mechanism.2  The  results  computed  with  this 
corrected  mechanism  were  essentially  the  same  as  those  obtained  in  the  preliminary  runs:4  the 
formation  of  polycyclic  aromatic  hydrocarbons  (modeled  in  our  work  up  to  the  size  of  coronene, 
C24H12)  via  the  ionic  reaction  pathway  is  much  slower  than  that  via  the  pathway  involving 
neutrals.  The  factors  responsible  for  that  are:4  (1)  a  relatively  low  rate  of  chemi-ionization  as 
compared  to  that  of  radical  production  at  the  same  conditions,  and  (2)  reversibility  of  the  reaction 
steps  responsible  for  the  formation  of  high-molecular  weight  species.  The  latter  condition  is  what 
make  the  reaction  network  so  complex,  to  the  extent  that  no  simple  analysis  based  on  equilibrium 
assumptions  or  irreversible  kinetics  can  be  applied. 

We  also  performed  flame  simulations  with  a  newer  ion  reaction  mechanism  of  AeroChem.  The 
results,  as  far  as  the  ions-vs-neutrals  competition  is  concerned,  were  again  the  same  as  those 
obtained  with  the  previous  mechanisms,  and  discussed  above.  These  results,  however,  cannot  be 
taken  as  final,  because  the  ion  reaction  mechanism  is  not  yet  completely  updated  and  tested.1  For 
example,  in  the  present  state,  it  produces  an  order  of  magnitude  higher  total  concentration  of  ions 
than  determined  in  experiment.5  This  point  and  others  will  have  to  be  resolved  in  the  future  work. 

Neutral  Reaction  Mechanism 

Accurate  prediction  of  polycyclic  aromatic  hydrocarbons  (PAHs)  in  flames  is  critical  for 
modeling  of  soot  formation  in  hydrocarbon  combustion.  In  our  initial  efforts,  the  objective  was  a 
proof-of-a-concept  computer  simulation  of  PAH  production  in  high-tempenture  hydrocarbon 
pyrolysis,3  and  later,  in  a  low-pressure  laminar  premixed  acetylene  flame.6  Now,  when  the 
acetylene-addition  mechanism  we  have  advocated  for  the  last  five  years  is  well  tested  under  various 
experimental  conditions  and  by  different  research  groups,  the  objective  is  to  improve  the 
mechanism  to  a  point  that  it  is  capable  of  quantitative  predictions.  Preliminary  results  of  such 
efforts  are  reported  below. 
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We  chose  to  model  the  PAH  profiles  in  three  well-studied  laminar  premixed  flames,  whose 
conditions  are  summarized  in  the  following  table: 


Table  1.  Flame  conditions 


No. 

Composition  in  Argon 

P  (torr) 

Cold  gas  velocity  (cm/s) 

Reference 

1 

46.5%  C2H2-48.5%C)2 

50 

MIT  group7 

2 

23.6%  C2H2  —  21.4%  O2 

90 

20 

Bockhom  et  al.  8 

3 

16.5%  C2H4  -  17.9%  O2 

760 

7.8 

Harris  et  al? 

These  flames  were  computed  using  the  Sandia  burner  code.10  The  Lennard- Jones  parameters  of 
PAH  molecules  were  calculated  using  empirical  correlations  and  available  physical  properties.11 
The  transport  coefficients  obtained  in  this  manner  agreed  well  with  experimental  data  available  for 
several  compounds.  The  thermodynamic  data  were  taken  from  Stein  and  Fahr,12  except  for  the 
Cf-H  group  additivity  which  was  based  on  the  recent  results  of  Green  et  a/.13 

The  reaction  mechanism  is  composed  of  two  parts:  the  small  molecule  reactions,  responsible 
for  the  main  flame  structures;  and  reactions  describing  the  formation  and  growth  of  PAHs.  The 
reaction  mechanism  for  the  former  part  was  composed  of  data  taken  primarily  from  several 
sources:  (a)  our  recently  optimized  methane  oxidation  mechanism,14  (b)  revised  mechanism  of 
Frenklach  and'Wamatz,6  and  (c)  revised  acetylene  oxidation  mechanism  of  Hwang  et  a/.15  The 
main  revisions  included:  the  rate  coefficients  of  acetylene-hydroxyl,  and  acetylene-oxygen  atom 
reactions 

C2H2  +  OH  -»  products 
C2H2  +  O  -» products 

were  taken  from  Miller  and  Melius,1^  and  Mahmud  and  Fontijn,17  respectively.  The  rate 
coefficient  for  reaction 

C2H2  +  C2H3  ->  C4H4  +  H 
was  taken  from  Fahr  and  Stein.1* 

The  formation  of  the  first  aromatic  ring  followed  the  basic  reaction  scheme  of  Frenklach  and 
Wamatz,6 

n-C4H3  +  C2H2  — »  phenyl 
n-QHs  +  C2H2  -»  benzene  +  H 
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Mole  fraction 


however  with  the  rate  coefficients  taken  from  the  work  of  Westmoreland  et  a/.,19  who  used  a 
QRRK  methodology  of  estimating  chemical  activation  processes  involved  in  these  reactions.  The 
assignment  of  the  rate  coefficients  for  the  main  reaction  classes  responsible  for  PAH  growth  are 
summarized  in  Table  2. 


Table  2.  Revision  of  rate  coefficients  for  main  reaction  classes  of  PAH  growth 


Reaction 

Prototype  reaction 

Reference 

aryl-H  +  H  -»  aryl  +  H2 

benzene  +  H  — » phenyl  +  H2 

Kiefer  et  al -20 

aryl  +  C2H2  — >  aiyl-C2H  +  H 

phenyl  +  C2H2  — » phenyl-acetylene  +  H 

Fahr  and  Stein18 

aryl  +  O2  — >  products 

phenyl  +  O2  ->  products 

Lin  and  Lin21 

The  computational  results  for  the  three  flames  are  presented  in  Figures  1  to  6.  These  results 
were  obtained  with  the  rate  coefficients  for  several  small-molecule  reactions  adjusted  by  a  factor  of 
2  to  3,  still  within  their  experimental  uncertainties.  As  can  be  seen  in  these  figures,  the  agreement 
between  the  modeling  predictions  and  experimental  data  is  good  for  major  species.  The  agreement 
for  PAHs  is  not  exactly  quantitative,  but  certainly  much  closer  —  in  both  shape  and  absolute 
values  —  than  reported  previously.6  The  important  result  here  is  that  the  near-quantitative 
accuracy  was  obtained  for  PAHs  with  a  single  reaction  mechanism  for  all  the  three  flames 
simultaneously.  The  major  reaction  pathways,  their  roles  in  PAH  growth,  and  the  sensitivity 
information  are  essentially  the  same  as  reported  and  discussed  in  the  previous  work.6 
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Figure  1.  Computed  vs.  experimental  major 
species  concentrations  of  Flame  1. 


Figure  2.  Computed  vs.  experimental  major 
aromatic  species  concentrations  of  Flame  1 . 
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Figure  3.  Measured  mole  fractions  of  diacetylene 
and  selected  PAHs  of  Flame  2. 
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Figure  4.  Computed  mole  fractions  of 
diacetylene  and  PAHs  of  Flame  2. 


Figure  5.  Computed  vs.  experimental  major 
species  concentrations  of  Flame  3. 


Figure  6.  Computed  vs.  experimental 
minor  species  concentration  of  Flame  3. 


Further  work  on  improving  the  prediction  of  PAH  profiles  under  the  flame  conditions  and 
analysis  of  the  factors  affecting  it  is  in  progress. 
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Using  the  neutral  reaction  mechanism  described  in  the  previous  section,  we  performed  detailed 
modeling  of  soot  particle  nucleation  and  growth  in  several  laminar  premixed  flames.  The 
computational  model  consists  of  three  logical  parts:  (I)  initial  PAH  formation ,  which  includes  a 
detailed  chemical  kinetic  description  of  acetylene  pyrolysis  and  oxidation,  formation  of  the  first 
aromatic  ring,  and  its  subsequent  growth  to  a  prescribed  size;  (II)  planar  PAH  growth,  comprised 
of  replicating-type  growth  of  PAHs  beyond  the  prescribed  size;  and  (III)  spherical  particle 
formation  and  growth,  consisting  of  coagulation  of  PAHs  formed  in  part  (II)  followed  by  the 
growth  of  the  resulting  particles  by  coagulation  and  surface  reactions. 

The  initial  part  of  the  particle  dynamics  model,  nucleation,  describes  the  planar  growth  of 
PAHs  via  the  H-abstraction/C2H2-addition  reaction  sequence3-6  using  the  technique  of  chemical 
lumping.22*23  This  method  provides  a  mathematically  rigorous  description  of  the  growth  process 
to  an  infinite  size  PAH,  which  can  be  schematically  represented  as 

A  /  -4  A /+;  -4  A/+2  -4  —  -4  Aoo  ,  (1) 

where  A/  represents  an  aromatic  species  containing  /  fused  rings.  The  PAH  species  formed  in  (1) 
are  then  allowed  to  coagulate,  that  is,  all  the  A,-’s  (j  =  /,  1+1,...,°°)  collide  with  each  other  forming 
dimers;  the  dimers,  in  turn,  collide  with  A/  forming  trimers  or  with  other  dimers  forming  tetramers; 
and  so  on.  The  coalescence  reactions  were  treated  as  irreversible  having  sticking  coefficients  of 
unity.  As  the  focus  of  this  work  was  on  very  young,  small  particles,  it  was  assumed  that  the 
coagulation  dynamics  is  in  the  free-molecular  regime.  A  size-independent  enhancement  factor  of 
2.2,  based  on  the  results  of  Harris  and  Kennedy,24  was  used  in  calculations  of  collision 
frequencies. 

Beginning  with  the  dimers,  the  forming  clusters  were  assumed  to  be  “solid  phase”  and  allowed 
to  add  and  lose  mass  by  surface  reactions 


SOOt- 

■H 

+  H  <=>  Csoot-  +  H2 

(2) 

SOOt* 

+ 

H  -4  Cgoof-H 

(3) 

soot* 

+ 

C2H2  -4  C^f-H  +  H 

(4) 

soot* 

+ 

O2  “4  products 

(5) 

soot- 

H 

+  OH  -4  products , 

(6) 
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where  Csoot-H  represents  an  arm-chair  site  on  the  soot  particle  surface  and  Csoot 
corresponding  radical.  This  mechanism  is  adopted  based  on  the  postulate25'27  that  the  H- 
abstraction/C2H2*addition  (HACA)  reaction  sequence3  is  responsible  for  high-temperature  growth 
of  all  forms  of  carbonaceous  materials.  Wieschnowsky  et  al.2%  also  suggested  that  the  HACA 
reactions3*6  “offer  a  key  to  the  understanding  of  a  number  of  phenomena”  observed  in  their  flame 
study. 

Estimation  of  the  rate  coefficients  of  the  heterogeneous  reactions  was  based  on  analogous  gas- 
phase  reactions  of  one-ring  aromatics,  benzene  and  phenyl.  The  reference  gaseous  reactions  for 
heterogeneous  reactions  (2>— (5)  were  the  corresponding  gas-phase  reactions  in  Table  (2). 
Oxidation  by  OH,  reaction  (6),  was  assumed  to  proceed  with  the  collision  efficiency  of  0.13, 
based  on  the  results  of  Neoh  et  al 29  The  particle  dynamics  —  the  evolution  of  soot  particles 
undergoing  simultaneous  nucleation,  coagulation  and  surface  reactions  described  above  —  was 
modeled  by  a  method  of  moments  (Method  II  of  Ref.  30).  This  method  does  not  require  the 
assumption  of  a  particle  size  distribution  function  (PSDF).  The  closure  of  the  differential 
equations  for  the  PSDF  moments  is  accomplished  by  interpolation  between  the  moments.22*30  The 
numerical  integration  of  the  moment  equations  took  from  20  to  50  s  on  an  IBM  3090/600S  main¬ 
frame  computer. 

Figures  7  and  8  present  the  results  computed  using  the  model  discussed  above  for  two  of  the 
chosen  flames:  Flame  B  —  flame  No.  1  of  Bockhom  and  co-workers28  —  25.4%  C2H2  -  19.6 
%  O2  -  Argon,  pressure  90  torr;  and  Flame  H  —  an  atmospheric  16.5%  C2H4  -  17.9%  O2  - 
Argon,  of  Harris  and  co-workers. 31  Both  flames  were  computed  with  the  same  kinetic  model, 
assuming  A4  (pyrene  or  acepyrene)  to  be  the  A /  species  which  initiates  the  planar  growth  of  PAHs 
and  the  formation  of  particles,  parts  (II)  and  (III)  of  our  model,  respectively.  Changing  A /  affects 
strongly  the  rate  of  particle  nucleation. 

It  can  be  seen  in  Figs.  7  and  8  that  the  model  predictions  are  in  relatively  close  agreement  with 
experiment  for  the  initial,  particle  inception  part  of  the  flames;  however,  the  computed  and 
experimental  results  begin  to  deviate  from  one  another  with  increasing  flame  height.  One  factor 
identified  as  a  cause  of  this  disagreement  is  the  shape  of  the  PAH  profiles.  This  is  demonstrated  in 
Fig.  7,  where  it  is  shown  that  using  a  narrower  (experimental-like8)  profile  for  A4  improves  the 
agreement  between  the  model  and  experiment  for  Flame  B.  It  was  also  computed  that  PAH 
condensation  on  the  particle  surface  is  sufficiently  fast  to  significantly  affect  the  PAH  profiles  in 
the  post-flame  zone. 
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Fig.  7.  Comparison  of  model  predictions  and  experimental  data  for  Flame  B:  (a)  particle  number  density;  (b) 
particle  surface  area;  (c)  average  particle  diameter  [circles  -  experimental  results,^®  solid  lines  -  model  prediction, 
dotted  lines  -  model  prediction  using  an  experimental-like  profile  for  A4];  (d)  specific  rale  of  soot  particle  oxidation 
by  O2  [solid  line  -  model  prediction,  dashed  line  -  from  the  expression  of  Nagle  and  Strickland-Constable.33 
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Fig.  8.  Comparison  of  model  predictions  and  experimental  data  for  Flame  H:  (aMc)  N  -  particle  number  density,  r 
-  particle  radius,  <rp>  -  average  d  -  average  particle  diameter,  circles  -  experimental  results,31  solid  lines  - 
model  prediction;  (d)  specific  surface  growth  rate  coefficient. 
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As  expected  from  the  dynamics  of  Smoluchowski  coagulation,  the  computed  rate  of  nucleation 
is  balanced  by  the  rate  of  coagulation  throughout  the  particle  inception  zone.  Particle  inception  is 
primarily  determined  by  PAH  coagulation,  initiated  and  controlled  by  PAH  coalescence  into 
dimers.  Our  model  describes  the  surface  processes  in  terms  of  elementary  reactions  of  active  sites. 
The  density  of  active  sites  is  determined  by  the  gas-phase  environment  —  by  the  competition 
among  reactions  (2)  -  (6).  Hence,  the  dissimilarity  in  surface  growth  behavior  can  be  explained  by 
the  difference  in  the  dynamics  of  active  sites.  Our  model  resolves  the  controversy  brought  up 
recently  by  Bockhom  and  co-workers,28  who  observed  that  the  surface  growth  rate  is  independent 
of  the  surface  area  of  soot  particles,  contrary  to  Harris  and  Weiner,32  who  postulated  that  the  rate 
of  surface  growth  is  first  order  in  particle  surface  area.  Also,  the  computed  net  surface  growth  rate 
is  computed  to  be  in  close  agreement  with  that  determined  by  Harris  and  Weiner32  (Fig.  8d)  and 
the  predicted  rate  of  soot  oxidation  by  O2  agrees  well  with  the  expression  of  Nagle  and  Strickland- 
Constable34  (Fig.  7d). 

Our  model  predicts  what  is  essentially  the  classical  picture  of  soot  particle  inception:34 
formation  of  PAHs,  their  coagulation  into  clusters  and  cluster  coagulation,  and  surface  growth  due 
to  chemical  and  physical  processes.  The  model  also  predicts  the  classical  structure  of  soot 
particles:35  a  less  dense  particle  core,  composed  of  randomly  oriented  PAH  oligomers,  and  a  more 
dense  concentrically-arranged  particle  shell. 


FUTURE  WORK 

The  objectives  for  the  next  year  are  the  following: 

1 .  To  complete  the  revision  and  testing  of  the  neutral  reaction  mechanism  of  PAH  production,  and 
to  prepare  a  manuscript  reporting  these  results. 

2.  To  apply  the  developed  mechanism  and  computer  code  for  soot  particle  nucleation  and  growth 
to  other  combustion  environments  of  interest. 

3.  To  review  the  ionic  reaction  mechanisms,  which  will  be  supply  to  us  by  AeroChem. 

4.  To  perform  flame  simulations  with  the  ionic  mechanism(s)  of  (3)  above. 

5 .  To  analyze  the  computational  results  of  (4)  above. 
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PRESENTATIONS  AND  PUBLICATIONS 

1 .  Frenklach,  M.,  “Particulate  Formation  in  Hydrocarbon  Combustion.’’  Invited  seminar  at  the 
University  of  Cincinnati,  Department  of  Chemical  and  Nuclear  Engineering,  February  2, 1989. 

2.  Frenklach,  M.,  “Formation  of  PCA,”  Prepared  statement  at  the  workshop  on  “Current 
Problems  in  Soot  Formation  during  Combustion,  Especially  the  Mechanism  of  Soot 
Formation,”  the  Commission  on  Condensation,  Academy  of  Sciences  Gottingen,  FRG,  March 
29-30,  1989. 

3.  Frenklach,  M.,  “Surface  Growth,”  Prepared  comment  at  the  workshop  on  “Current  Problems 
in  Soot  Formation  during  Combustion,  Especially  the  Mechanism  of  Soot  Formation,”  the 
Commission  on  Condensation,  Academy  of  Sciences  Gottingen,  FRG,  March  29-30, 1989. 

4.  Wang,  H.  and  Frenklach,  M.,  “Modeling  of  PAH  Profiles  in  Premixed  Flames”.  Paper 
presented  at  the  Fall  Technical  Meeting  of  the  Eastern  States  Section  of  the  Combustion 
Institute,  Albany,  New  York,  October  1989,  Paper  12. 

5.  Frenklach,  M.  and  Wang,  H.,  “Detailed  Modeling  of  Soot  Particle  Nucleation  and  Growth.” 
Paper  submitted  to  the  Twenty-Third  International  Symposium  on  Combustion. 

6.  Frenklach,  M.,  “Reduction  of  Chemical  Reaction  Models.”  In  Numerical  Approaches  to 
Combustion  Modeling  (E.  S.  Oran  and  J.  P.  Boris,  Eds.),  Am.  Inst.  Aeronautics  Astronautics 
Prog.  Ser.,  in  press. 


PROFESSIONAL  PERSONNEL 

Dr.  Michael  Frenklach  —  Principal  investigator 
Mr.  Hai  Wang  —  Graduate  student 


None. 
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